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Abstract 

We first consider the Hamiltonian formulation of n = 3 systems in 
general and show that all dynamical systems in R'^ are bi-Hamiltonian. 
An algorithm is introduced to obtain Poisson structures of a given 
dynamical system. We find the Poisson structures of a dynamical 
system recently given by Bender et al. Secondly, we show that all 
dynamical systems in M" are {n — l)-Hamiltonian. We give also an 
algorithm, similar to the case in M^, to construct a rank two Poisson 
structure of dynamical systems in M". We give a classification of the 
dynamical systems with respect to the invariant functions of the vector 
field X and show that all autonomous dynamical systems in M" are 
super-integrable. 







1. Introduction 



Hamiltonian formulation of n = 3 systems has been intensively considered 
in the last two decades. Works [1], [2] on this subject give a very large 
class of solutions of the Jacobi equation for the Poisson matrix J. Recently 
generalizing the solutions given in [1] we gave the most general solution of 
the Jacobi equation in R^, [3]. Matrix J = (J*-'), i,j = 1, 2, • • • , n defines 
a Poisson structure in M" if it is skew-symmetric, J'^^ — — J^*, and its entries 
satisfy the Jacobi equation 

f'di J^" + f^di J^' + f^di = 0, (1) 

where i, j, A; = 1, 2, ■ ■ ■ , n. Here we use the summation convention, meaning 
that repeated indices are summed up. We showed in [3] that the general 
solution of the above equation (1) in the case n = 3 has the form 

f^^^x^^^^k^, = 1,2,3, (2) 

where and ^ are arbitrary differentiable functions of i — 1,2,3 and 
e^^^ is the Levi-Civita symbol. Here t should be considered as a parameter. 
In the same work we have also considered a bi- Hamiltonian representation of 
Hamiltonian systems. It turned out that any Hamiltonian system in has 
a bi-Hamiltonian representation. 

In the present paper we prove that any n-dimensional dynamical system 

f (3) 

where ), is Hamiltonian, that is, has the form 

x^^f^djH, i = l,2,...,n, (4) 

where J = (T^) is a Poisson matrix and H, as well as J'-' , are differentiable 
functions of the variables xi,X2, ■ ■ ■ , Xn, t. Moreover, we show that the system 
(3) is (n — l)-Hamiltonian. This problem in the case n = 3 was considered 
in [4], [5] where authors start with an invariant of the dynamical system as 
a Hamiltonian and then proceed by writing the system in the form (4) and 
imposing conditions on J so that it satisfies the Jacobi equation. But proofs 
given in these works are, as it seems to us, incomplete and not satisfactory. 
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Using (2) for matrix J we can write equation (4) in as 



f=/iV*xVi/. (5) 

Let X be a vector field in M?. If Hi and H2 are two invariant functions of 
X, i.e., X{H^) = X^djH^ = 0, a = 1, 2, tlien X is parallel to Vifi x VH2. 
Therefore 

X = nVHi X V//2, (6) 

where the function is a coefficient of proportionality. The right-hand side 
of equation (6) is in the same form as the right-hand side of equation (5), so 

— * 

X is a Hamiltonian vector field. We note that the equation which allows to 
find the invariants of a vector field X is a first order linear partial differential 
equation. We remark here that dynamical systems in differ from the 
dynamical systems in for n > 3. We know the general solution (2) of the 
Jacobi equation (1) in W. In R", as we shall see in the last section, we know 
only the rank 2 solutions of the Jacobi equations for all n. 

An important difference of our work, contrary to other works in the sub- 
ject, is that in the construction of the Poisson structures we take into account 
the invariant functions of the vector field X rather than the invariants (con- 
stants of motion) of the dynamical system. The total time derivative of a 
differentiable function F in along the phase trajectory is given by 

'-f = f + X.VF. (7) 

dt dt ^ ^ 

— * — * — * 

An invariant function of the vector field X{x , x , . . . ,x",t), i.e., X • VF = 0, 
is not necessarily an invariant function (constant of motion) of the dynamical 
system. For autonomous systems where X = X( ) these invari- 

ant functions are the same. We give a representation of the vector field X 
in terms of its invariant functions. We show that all autonomous dynamical 
systems are super-integrable. A key role plays the existence of n — 1 func- 
tionally independent solutions C,a{x^ i ^'^ i ■ ■ ■ , 0; (o; = 1, 2, • • • , n — 1) of the 
linear partial differential equation 

X.VC.X'§ + X^§ + ... + X"^ = 0, (8) 
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where X* = X*(a;^, a;^, . . . , t), i = 1,2, ■ ■ ■ ,n, are given functions (see 
[6]- [8]). For all a = 1, 2, • ■ ■ , n — 1, VCa is perpendicular to the vector field 

— * 

X. This leads to the construction of the rank 2 Poisson tensors for n > 3: 

Z*-? = //f'^°l"2""n-2 ,Wl-in-2 a. R C ■■■8 C fQl 

<^a c ^jiS,ai Sa2 K7n-2S,a„_2) \'^J 

where j = 1, 2, • • • , n, and a = 1, 2, • • • , Here e^^'^i-^"-^ and e"«i«2 - an-2 
are Levi-Civita symbols in n and n — 1 dimensions respectively. Any dynam- 
ical system with the vector field X possesses Poisson structures in the form 
given in (9). Hence we can give a classification of dynamical systems in 

— * 

M" with respect to the invariant functions of the vector field X. There are 
mainly three classes where the super-integrable dynamical systems consti- 
tute the first class. By the use of the invariant functions of the vector field 
X{x^,x'^, . . . ,x"',t) in general we give a Poisson structure in which has 
rank 2. For autonomous systems, the form (9) of the above Poisson structure 
first was given in the works [11] and [12]. 

Our results in this work are mainly local. This means that our results 
are valid in an open domain of M" where the Poisson structures are different 
from zero. In [3] we showed that the Poisson structure (2) in preserves 
its form in the neighborhood of irregular points, lines and planes. 

In the next section we give new proofs of the formula (2) and prove 
that any dynamical system in is Hamiltonian. So, following [3] we show 
that any dynamical system in is bi-Hamiltonian. Applications of these 
theorems to several dynamical systems are presented. Here we also show 
that the dynamical system given by Bender at al [10] is bi-Hamiltonian. In 
section 3 we discuss Poisson structures in M". We give a representation of 
the Poisson structure in R" in terms of the invariant functions of the vector 

— * 

field X. Such a representation leads to a classification of dynamical systems 
with respect to these functions. 

2. Dynamical Systems in 

Although the proof of (2) was given in [3], here we shall give two simpler 
proofs. The first one is a shorter proof than the one given in [3]. In the 
sequel we use the notations — x,x^ — y, x^ — z. 

Theorem 1. All Poisson structures in R^ have the form (2), i.e., J^^ — 
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fie'^'^ dk Hq. Here fi and Hq are some differentiable functions of and t, 
Cz = 1,2,3; 

Proof. Any skew-symmetric second rank tensors in can be given as 

f^^e'^'^Jk, = 1,2,3, (10) 

where Ji , J2 and J3 are differentiable functions in and we assume that 
there exists a domain in so that these functions do not vanish simulta- 
neously. When (10) inserted into the Jacobi equation (1) we get 

J. (Vxj)-0, (11) 

where J = (Ji, J2, J3) is a differentiable vector field in R"^ not vanishing in 
fl. We call J as the Poisson vector field. It is easy to show that (11) has 

— * — * 

a local scale invariance. Let J — ^ E, where ip is an arbitrary function. If 
E satisfies (11) then J satisfies the same equation. Hence it is enough to 
show that E is proportional to the gradient of a function. Using freedom of 
local scale invariance we can take E = {u,v, 1) where u and v are arbitrary 
functions in R^. Then (11) for vector E reduces to 

dyU — dxV — vdzU -\- udzV — 0, (12) 

where x,y,z are local coordinates. Letting u — ^ and v = where / 
and p are functions oi x,y, z we get 

dj dy{p - dj) - djd,{p - dj) = 0. (13) 

General solution of this equation is given by 

p-dj^h{f,z), (14) 

— * 

where h is an arbitrary function of / and z. Then the vector filed E takes 
the form 

E = {dJ, dyf, dJ + h). (15) 

Let g{f, z) be a function satisfying Q z = hdfQ. Here we note that g{f, z) 

dg 



% dzf + g,z where g^^ = ds g{f{x, y, z), s)\s=z- Then (15) becomes 



E- .r,r^,,r, Vg, (16) 
{dJ + h)dfg 
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which completes the proof. Here df g = □ 

The second proof is an indirect one which is given in [8] (Theorem 5 in this 
reference) . 

Definition 2. Let F be a vector field in M . Then the equation F ■ dx — 

is called a Pfaffian differential equation. A Pfaffian differential equation is 
called integrahle if the 1-form F ■ dx = /idH, where /i and H are some 
differentiable functions in . 

Let us now consider the Pfaffian differential equation with the Poisson vector 
field J in (10) 

J-dx = 0. (17) 
For such Pfaffian differential equations we have the following result (see [8]). 

Theorem 3. A necessary and sufficient condition that the Pfaffian differen- 
tial equation J ■ dx = should be integrable is that J • (V x J) = 0. 

— * — » 

By (11), this theorem implies that J = /iV\l/. 

A well known example of a dynamical system with Hamiltonian structure 
of the form (4) is the Euler equations. 

Example 1. The Euler equations [6] are 

h - h 
X = ———yz, 

y = ^ J ^ xz, (18) 

^1 — -12 

where Ii,l2,h £ R are some (non- vanishing) real constants. This system 
admits Hamiltonian representation of the form (4). The matrix J can be 
defined in terms of functions ^ — Hq — —\{x'^ + + z^) and — 1, and we 

v'^ z^ 

Writing the Poisson structure in the form (2) allows us to construct bi- 
Hamiltonian representations of a given Hamiltonian system. 



5 



Definition 4. Two Poisson structures Jq and Ji are compatible, if the sum 
Jo + Ji defines also a Poisson structure. 

Lemma 5. Let ^^Hq, and Hi he arbitrary differentiahle functions. Then the 
Poisson structures Jq and Ji given by Jq = ^e'^^^dk Hq and Ji = —fie^^^dk Hi 
are compatible. 

This suggests that all Poisson structures in have compatible companions. 
Such compatible Poisson structures can be used to construct bi-Hamiltonian 
systems (for Hamiltonian and bi-Hamiltonian systems see [6], [9] and the ref- 
erences therein). 

Definition 6. A Hamiltonian equation is said to he bi-Hamiltonian if it ad- 
mits compatible Poisson structures Jo and J\ with the corresponding Hamil- 
tonian functions Hi and Hq respectively, such that 

dx 

--^JoVHi = JiVHo. (19) 



Lemma 7. Let Jq be given by (2), i.e., .Jq = /xe'-'^Sfc Jfo, o-nd let Hi be any 
differentiahle function, then the Hamiltonian equation 

— = JoVHi = Ji V Ho ^11 WHi X Vi^o, (20) 

(JjV 

is bi-Hamiltonian with the second Poisson structure given by Ji with entries 
J I = —fxe'^^'^dkHi and the second Hamiltonian Hq. 

Let us prove that any dynamical system in has Hamiltonian form. 

Theorem 8. All dynamical systems in are Hamiltonian. This means 
that any vector field X in M'^ is Hamiltonian vector field. Furthermore all 
dynamical systems in are bi-Hamiltonian. 

Proof. Let ( be an invariant function of the vector field X, i.e., X{() = 
X ■ = 0. This gives a first order linear differential equation in for ^. 
For a given vector field X = (/, g, h) this equation becomes 

dC dC dC 

f{x, y, ^' ^) ^ + 9{x, y, z,t)-^ + h{x, y, z,t)-^ = 0, (21) 
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where x, y, z are local coordinates. From the theory of first order linear 
partial differential equations [6], [7], [8] the general solution of this partial 
differential equation can be determined from the following set of equations 

dx dy dz 

^ (22) 



f{.x,y,z,t) g{x,y,z,t) h{x,y,z,t)' 
There exist two functionally independent solutions Ci and C2 of (22) in an 
open domain D C M.^ and the general solution of (21) will be an arbitrary 
function of Ci and C2, i-e., C — ^{Ci: Ci)- This implies that the vector field X 
wiU be orthogonal to both VCi and VC2- Then X — /i (VCi) X (VC2)- Hence 
the vector field X is Hamiltonian by (5). □ 

This theorem gives also an algorithm to find the Poisson structures or the 
functions Hq, Hi and of a given dynamical system. The functions Hq and 
Hi are the invariant functions of the vector field X which can be determined 
by solving the system equations (22) and fi is determined from 

= /"^^ . (23) 

X-{VHoxVHi) 

Note that can also be determined from 

X^ 

^ d2HQd^Hi — dzHQd2Hi 

^2 



d^Hid^Hi - diHod^Hi 

)e 

diHod2Hi - d2HodiHi 



Example 2. As an application of the method described above we consider 
Kermac-Mckendric system 

X — —rxy, 

y = rxy -ay, (25) 
i = ay, 

where r, a G M are constants. Let us put the system into Hamiltonian form. 
For the Kermac-Mckendric system, equations (22) become 

dx dy dz 



-rxy rxy — ay ay 



(26) 
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Here a and r may depend on t in general. Adding the numerators and 
denominators of (26) we get 

dx dx + dy + dz , 

— . (27) 



—rxy 

Hence Hi = x + y + z is one of the invariant functions of the vector field 
Using the first and last terms in (26) we get 

dx dz 



-rx a 



(28) 



which gives Hq = r z+a In x as the second invariant function of the vector field 
X. Using (23) we get = xy. Since X — ixVHq x Vi^i, the system admits 
a Hamiltonian representation where the Poisson structure J is given by (2) 
with ^ — xy, — Hq — rz -\- a\n.x, and the Hamiltonian is Hi — x -\- y -\- z. 

Example 3. The dynamical system is given by 

X = yz{l + 2x'^ N/D), 

y = -2xz{l-y'^N/D), (29) 
z = xy{l + 2z^N/D), 

where N — x"^ + y^ + z^ — 1, D — x^y^ + y^z^ + A-X^z^. This example was 
obtained by Bender et all [10] by complexifying the Euler system in Example 
1. They claim that this system is not Hamiltonian apparently bearing in mind 
the more classical definition of a Hamiltonian system. Using the Definition 
6 we show that this system is not only Hamiltonian but also bi-Hamiltonian. 
We obtain that 



H,= ^^i^^N, Hi = ^^{2y^z' + Ax\' + y' + 2xY-y'). (30) 
Here 

4[3D2 + DP + g]' ^^^^ 

where 

P = -2a;^ + V - 4xV + - V - + 14^^ + '^^ 

Q = -2x^ + I2x^z^ + 2x^ - 20x^z^ - 6x^z^ - b2x^z^ - Qx^z^ 

+y^ -y^ + Ay^z^ - IGy^z^ - 2z^ + 2z^. (32) 
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Indeed these invariant functions were given in [10] as functions A and B. 
The reason why Bender et al [10] concluded that the system in Example 3 is 

— * 

non-Hamiltonian is that the vector filed X has nonzero divergence. It follows 
from X — iiVHq x V-ffi that V ■ = 0. When ji is not a constant the 

corresponding Hamiltonian vector field has a nonzero divergence. 

Remark 1. With respect to the time dependency of invariant functions of 
the vector field X dynamical systems in can be spht into three classes. 

Class A. Both invariant functions Hq and Hi of the vector field X do not 
depend on time explicitly. In this case both Hq and are also invariant 
functions of the dynamical systems. Hence the system is super-integrable. 
All autonomous dynamical systems such as the Euler equation (Example 1) 
and the Kermac-Mckendric system (Example 2) belong to this class. 

— * 

Class B. One of the invariant functions Hq and Hi of the vector field X 

depends on t explicitly. Hence the other one is an invariant function also of 
the dynamical system. When Ji, I2 and in Example 1 are time dependent 
the Euler system becomes the member of this class. In this case Hq is the 
Hamiltonian function and Hi is the function defining the Poisson structure. 
Similarly, in Example 2 we may consider the parameters a and r as time 
dependent. Then Kermac-Mckendric system becomes also a member of this 
class. 

Class C. Both Hq and Hi are explicit functions of time variable t but 
they are not the invariants of the system. There may be invariants of the 
dynamical system. Let F be such an invariant. Then 

where for any F and G 

{F, G}, = di FdjG, a = 0, 1 . (34) 



3. Poisson structures in R 

Let us consider the dynamical system 
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— ^X\x\x\---,x\t), z = l,2,---,n. 



(35) 



Theorem 9. All dynamical systems in M" are Hamiltonian. Furthermore 
all dynamical systems in MJ^ are {n — 1) -Hamiltonian. 

Proof. Extending the proof of Theorem 8 to consider the hnear partial 

differential equation (8). There exist n — 1 functionally independent solutions 
Ha, {a = 1, 2, ■ ■ ■ , n — 1) of this equation (which are invariant functions of 
the vector field X) [6]-[8]. Since X is orthogonal to the vectors VHa, {a — 
1, 2, • • • , n — 1), we have 



X = 11 



ei 62 
diHi d2Hi 



diH, 



doH„ 



n-l '-'2^^n-l 



6n 

dnH, 



dr,H, 



n-'-'n— 1 



(36) 



where the function is a coefficient of proportionality and el is n-dimensional 
unit vector with the ith coordinate 1 and remaining coordinates 0. Therefore 

X' = i2e'^'^'-^"-^d,,Hid,,H2 ■ ■ ■ d,„_,Hn-i. (37) 

Hence all dynamical systems (35) have the Hamiltonian representation 



= J^djHa, i = l,2, (no sum on a) 



(38) 



with 



pj — ,.fr<^aia2-a„-2 MJl-Jn-2 a, TT f) U . . . f) ■ ff ( 'iQ) 

■-'q — Alt t <<7i-nai <<72 -"a2 '^Jn-2-"an-2 ) K'^^J 

where i, j = 1, 2, • • • , n, a = 1, 2, • • • , n - 1. Here e^^'^i-^"-^ and e«°i«2-an-2 
are Levi-Civita symbols in n and n — 1 dimensions respectively. The function 
/J, can be determined, for example, from 



X' 



d^H^ 



dnH^ 



(40) 



d2Hn-l 



■ dnHn-l 
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It can be seen that the matrix with the entries J^J given by (39) defines 
a Poisson structure in M" and since 



Ja-VHp^O, a,/? =1,2, 



1, 



(41) 



with (3 ^ a, the rank of the matrix equals 2 (for all a = 1, 2, • • • , n — 1). 
In (38) we can take any of Hi, H2, ■ ■ ■ , Hn~i as the Hamilton function and 
use the remaining H^s in (39). Wc observe that all dynamical systems (35) 
in have n — 1 number of different Poisson structures in the form given 
by (39). The same system may have a Poisson structure with a rank higher 
than two. The following example clarifies this point. 



Example 4. Let 

Xi = X4, X2 = X3, Xs = -X2, X4 = -Xi. 

Clearly this system admits a Poisson structure with rank four 



(42) 



J 



/ 1 \ 

10 
0-100 

\ -1 y 



H ^ -{xl + xl + xl + xl). 



(43) 



The invariant functions of the vector field X = (2:4, 2:3, —X2, —Xi) are 

1 



Hi 



- iXi -|- ^2 -|- -|- X^), 



2 (•*"2 + ^3)' 



H2 

Hs = X1X3-X2X4. 



(44) 

(45) 
(46) 



Then the above system has three different ways of representation with the 
second rank Poisson structures 



= fxe'^'''dkH,diH2, H = Hs, 



(47) 
(48) 
(49) 



where fi{xiX2 + XSX4) = 1. These Poisson structures are compatible not only 
pairwise but also triple- wise. This means that any linear combination of 
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these structures is also a Poisson structure. Let J = aiJi + q;2</2 + as^a then 
it is possible to show that 

^y,e'^kidkHidiH2, (50) 
where Hi and H2 are linear combinations of Hi, H2 and i^s, 

Hi = Hi- ^H2, H2 = - a2Hs if ^2 7^ 0, (51) 

a2 

Hi = aiHi - a^H2, H2 = H2 if = 0. (52) 

Definition 10. A dynamical system (35) in is called super-integrable if 
it has n — 1 functionally independent first integrals (constants of motion). 

Theorem 11. All autonomous dynamical systems inW^ are super-integrable. 

Proof. If the system (35) is autonomous, then the vector field X does not 
depend on t explicitly. Therefore each of the invariant functions H^, (a = 

— * 

l,2,---,n — l)of the vector field X is a constant of motion of the system 
(35). 

Some (or all) of the invariant functions H^, {a = 1,2, ■ ■ ■ ,n — 1) of the 
vector field X may depend on t. Like in we can classify the dynamical 
systems in M" with respect to the invariant functions of the vector field 
X{x^,x'^,---,x'',t). 

Class A. All invariant functions H^, (« = 1,2, ■■-.n — 1) of the vector 
field X do not depend on t explicitly. In this case all functions H^, (a = 
1, 2, • • • .n — 1) are also invariant functions (constants of motion) of the dy- 
namical system. Hence the system is super-integrable. In the context of the 
the multi- Hamiltonian structure, such systems were first studied by [12] and 
[11]. The form (39) of the Poisson structure was given in these works. Its 
properties were investigated in [13]. 

Clciss B. At least one of the invariant functions Ha, (a = 1, 2, ■ ■ • .n — 1) 
of the vector field X does not depend on t explicitly. That function is an 
invariant function also of the dynamical system. 
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Class C. All Hct^ (a = 1, 2, ■ ■ ■ .n — 1) are explicit function of time variable 
t but they are not the invariants of the system. There may be invariants of 
the dynamical system. Let F be such an invariant. Then 



dF ^ dF 
Itt ^~dt 
, where for any F and G 



(53) 



{F,G}a = J'JdiFdjG, a = 0,l,---,n-l. 



(54) 
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